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Abstract 


I x = I 0 e~^ x + B 


A technique for monitoring variability in the nu- 
clear absorption characteristics of aviation fuels has 
been developed. It is based on a highly collimated 
low energy gamma radiation source and a sodium 
iodide counter. The source and the counter assem- 
bly are separated by a geometrically well-defined test 
fuel cell. A computer program for determining the 
mass attenuation coefficient of the test fuel sample, 
based on the data acquired for a preset counting pe- 
riod, has been developed and tested on several types 
of aviation fuel. 

Introduction 

We have recently demonstrated (ref. 1) the fea- 
sibility of a nuclear gauging system for fuel quan- 
tity measurement onboard an aircraft. It is based 
on monitoring the number of photons arriving at 
a counting station located a predetermined distance 
from an Am 241 source capsule. Several such source 
capsule-counting station assemblies are judiciously 
located throughout the fuel tanks. When the en- 
tire pathlength between the source capsule and the 
counter station is occupied by the fuel, the number 
of surviving photons will be minimal. If, on the other 
hand, there is no fuel in the photon path, the num- 
ber of surviving photons will be maximum. By com- 
bining the information about the counting rates and 
the geometrical locations of the source-counter as- 
semblies in the tanks, a reliable measure of the total 
fuel content of the tank can be obtained at any time. 
This, of course, is true only if the quality of the fuel 
remains constant. Recently, the Airlines Electronic 
Engineering Committee (AEEC) has reported con- 
cern about the variability of aviation fuel character- 
istics as a function of the season and geographical 
origin (refs. 2, 3, and 4). The concern arises from 
variations in the fuel composition as well as the na- 
ture and amount of contaminants. In an effort to 
identify the degree of variability in fuel quality, we 
have set up a fuel characteristics monitoring system. 
This system and the computational procedure used 
to measure changes in the nuclear absorption coeffi- 
cients of the fuel samples are described in this report. 

Fuel Characteristics Monitoring System 

The fuel characteristics monitoring system is 
made up of a highly collimated 10 /iCi Am 241 
(59.5 keV) gamma radiation source and a 2-in- 
diameter x 2-in. Nal(Tl) crystal mounted on a photo- 
multiplier. The source and the counter assembly are 
separated by a 2-in-diameter x 4-in. glass fuel cell. 
The number of photons arriving at the Nal(Tl) crys- 
tal depends on the quality of the fuel in the fuel cell. 


wvhere 

F x number of photons arriving at the 

Nal (Tl) crystal 

f Q number of photons incident on the 

fuel cell 

t linear attenuation coefficient of the 

fuel 

k fuel path length (fuel cell length) 

8 background count (counts recorded 

in the absence of the source) 

By using a well-characterized medium in the fuel 
ell-— such as air or distilled water — the value of I Q 
an be determined from a measured value of I x . Once 
r o is determined for a fixed source-detector assem- 
bly, I x becomes the critical measurable parameter in 
he fuel quality study. An independent measurement 
■if the density of the test fuel, coupled with a value 
>f linear attenuation coefficient determined from the 
preceding equation, then permits a direct computa- 
ion of mass attenuation coefficient (/i/p, cm 2 /g) of 
he sample (where p is the density of the test fluid). 

Computational Procedure 

Basic description . A computer program for 
;amma radiation mass attenuation coefficient 
PGRMAC) has been written in MS-FORTRAN 77 
Version 3.31 for personal computers with a fixed-disk 
ystem, using MS-DOS Version 3.3. The program re- 
quires 12 156 bytes of disk space for storage. 

The program models the experimental procedure 
or calculating gamma ray attenuation coefficients 
n the test medium. The geometrical details of the 
est system are summarized in figure 1. Figure 2 is a 
>hotograph of the experimental system. As shown in 
tgure 1, gamma rays have to pass through air, glass 
uel cell walls, the test fluid, and a thin aluminum 
lousing to reach the detector surface. The intensity 
number of photons) of gamma radiation arriving at 
he detector can be written as follows: 

I x = I 0 {e~^ •^airl g /^gl ass ^g lass 

where 

x intensity of gamma radiation 

arriving at the detector 

o intensity of gamma radiation 

incident on fuel cell 



Mair 

Mglass 

linear attenuation coefficient for air 

linear attenuation coefficient for 
glass 

coefficient of the test fluid from equation (1). When 
the glass cell is filled with the test fluid, equation (1) 
can be written as 

Mai 

linear attenuation coefficient for 
aluminum 

Ix — Io ex P ( — Mair^airl ~ Mglass^glass 

— Mfluid-^fluid — Mal^al) (**) 

Mfluid 

linear attenuation coefficient for the 
test fluid 

^airl 

air path length 


x glass 

glass path length 

Then 

^al 

aluminum path length 


Mfluid 

B 

fluid path length 

background count (counts recorded 
in the absence of the source) 

Mfluid ~ [l n (7o/7 x ) — Mair^airl 

“ Mglass^glass — Mal^alj/^fluid (®) 


In equation (1), /4 a j r , ^airli Mglass> ^glass? Map ^ab 
and Xfl u jd are known parameters (refs. 1, 5, 6, and 
7). The intensities I Q and 7 X need to be calculated or 
measured. To determine 7 0 , we can choose air as the 
reference medium. Equation (1), after subtracting 
the background count, can be written as follows: 

I x — B = / 0 ^e _ ^ airXairl e~^g}ass x glass 

X g — ^fiuid^fluid g — ^al^al \ 


The calculated value from equation (6) is the lin- 
ear attenuation coefficient of the test fluid. Since the 
mass attenuation coefficient of the medium is of more 
fundamental importance than the linear attenuation 
coefficient, the density of the test fluid should be 
determined beforehand independently. In this pro- 
gram, the density of the test medium, measured ex- 
perimentally, is read as an input value, and the mass 
attenuation coefficient is simply equal to the linear 
attenuation coefficient divided by the density, i.e., 


7 , . . — Tip Mair^air 1 p Mglass^glass 
J x(air) “ 1 o\ c c 

^ g — £*air x fluid(air) g — / i al I al ^ ^2) 

This gives the following relation for I Q : 

Io = Ix (3-ir) exp | //air ^airl *b ^fluid(eiir) 

“b M gl ass ^ glass ~b Mal^al 1 (^) 


7 o — 7 x (air) exp (Mair*^air2 ~b Mglass 2 ' glass ~b Mal^al) 


( 4 ) 


where 

a 'air2 = ^airl “b a: fluid(air) 

Since 7 x (air) can be determined experimentally, 7 0 is 
readily calculated from equation (4). 

Once I 0 has been obtained from equation (4), the 
program proceeds to compute the linear attenuation 


(m) 


mass 


(m)i inear 

Density 


( 7 ) 


The region of the gamma ray spectrum in which 
we are interested is selected, and it is marked as the 
region of interest (ROI) in figure 3. This includes 
most of the area under the total capture peak in 
the spectrum. The energy distribution of the gamma 
source forms a peak in the ROI, and the area under 
the peak in the ROI is interpreted as 7 n for each 
medium in the program. The background counts 
within the ROI have to be determined and then sub- 
tracted from the full spectral counts in the ROI be- 
fore it can be used to determine 7 X or 7 0 . The ROI 
in this program is selected to be from channel 241 
to channel 421. The peak falls at the middle of the 
region, as seen in figure 3; this range is quite ade- 
quate to monitor and analyze the characteristics of 
each test fluid. 

Program input , output , and usage . The program 
input is made of two parts containing several param- 
eters and related file names. The first part of the 


2 



input consists of 10 essential parameters that are as- 
signed as constants in the program. They are listed 
below: 


MUAIR=0.0002132 
MUGLS=0. 46654 
MUAL=0.6639 
XAIR1=5.080 

XAIR2=XAIR1 + XFLUID(AIR) 


XGLS=0.680 

XAL=0.079 

XFUEL=10.062 

LBEG=241 

LEND=421 


In the list above, MUAIR, MUGLS, and MUAL 
are the linear attenuation coefficients of air, glass, 
and aluminum, respectively (refs. 5, 6, and 7); XAIR, 
XGLS, XAL, and XFUEL are the air, glass, alu- 
minum, and test medium pathlengths, respectively. 
The parameters LBEG and LEND correspond to the 
beginning and ending channel numbers for the ROI. 
The second part of the input is an eight-character 
string, which is read interactively from the keyboard, 
and which corresponds to two binary file names and 
one specific record in a data base. These two binary 
files, which are generated from the multichannel an- 
alyzer (MCA), refer to the gamma spectra through 
air and the test fluid taken on the same day. Each 
record in the data base has nine terms of informa- 
tion: density, temperature, type of medium, etc., in 
it. 

The program generates four formatted output 
files. Two of them are data files, converted from 
the binary, spectra through air and the test medium, 
that list channel numbers and counts in each of 
them as shown in tables 1(a) and (b). Plots of the 
data in tables 1(a) and (b) are shown in figures 4 
and 5, respectively. From the data file of the air 
spectrum, the program calculates the value of I Q and 
determines the location of the centroid of the peak 
channel, MAXAIR, for air on that day. The program 
then computes I x from the converted data file of the 
test fluid and determines the location of the peak 
centroid, MAXFLUID. 

MAXAIR and MAXFLUID should coincide and 
should remain constant as long as the electronic sys- 
tem gain remains unchanged. Thus a determination 
of E " LT D serves as a system calibration check. 

Once I 0 and I x have been determined by the above 
procedure, the linear attenuation coefficient of the 
test fluid (MUFLUID) can be obtained from equa- 
tion (6). Subsequently, the mass attenuation coeffi- 
cient (ATCFM) is readily calculable. The other two 
output files, the printout and monitor display, pro- 
vide a record of the desired information about the 
test medium. 

Program flowchart and listings . The program flow- 
chart and listings are summarized in appendices A 
and B, respectively. 


Test of the Sensitivity of the System 

To test the sensitivity of the monitoring sys- 
em, mass attenuation coefficient measurements were 
nade in common salt solutions in w'ater containing 
different amounts of salt. The experimental results, 
along with the corresponding calculated values, are 
summarized in table II and illustrated in figure 6. It 
s apparent from the data shown in this table that the 
measured and calculated values of the mass attenu- 
ition coefficients for different solutions agree within 
L0.5 percent. 

Applications 

The program PGRMAC has been applied to the 
measurement of mass attenuation coefficients for sev- 
eral types of aviation fuel. These fuels have been in- 
vestigated previously and can thus provide a good 
test of the validity of the computational procedure. 
Results for JP-4, JP-5, and Jet A are summarized in 
tables 111(a), (b), and (c), respectively. 

The results for the three types of fuel are further 
consolidated in table IV. 

Concluding Remarks 

A simple technique for monitoring nuclear absorp- 
tion characteristics of aviation fuels, based on low 
energy gamma ray attenuation in the test fuels, has 
been developed. It has been tested on three types 
of aviation fuels. It is noted that the mass attenua- 
tion coefficients of the three fuels are almost equal, 
even though their linear attenuation coefficients and 
densities are slightly different. It would therefore 
appear that a simultaneous measurement of linear 
and mass attenuation coefficients and densities is a 
highly informative procedure. It is further noted that 
the values of mass attenuation coefficients calculated 
by using the present procedure agree with the pre- 
viously reported values to within ±1 percent. It is 
therefore concluded that the procedures developed 
here are quite adequate for monitoring variability of 
>0.5 percent in fuel absorption characteristics as a 
function of the season and geographical points of ori- 
gin of the fuels. 

The international aviation consortium has agreed 
to provide us fuel samples from various parts of the 
world over the next 12 months. The linear and mass 
attenuation coefficients for 59.5 keV Am 241 gamma 
rays, as well as densities of the fuel samples, will be 
measured to assess the fuel composition variability. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
September 26, 1988 
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Table I. Cont nued 


(a) Concluded 


Channel 

number 

Counts 

Channel 

number 

Counts 

Clannel 
nu nber 

Counts 

Channel 

number 

Counts 

221 

10 

276 

123 

131 

1082 

386 

101 

222 

17 

277 

140 

132 

1067 

387 

97 

223 

13 

278 

165 

133 

1039 

388 

107 

224 

18 

279 

176 

134 

1032 

389 

111 

225 

19 

280 

214 

135 

1032 

390 

105 

226 

13 

281 

191 

136 

1017 

391 

76 

227 

1 

282 

227 

137 

1035 

392 

62 

228 

3 

283 

232 

138 

939 

393 

56 

229 

12 

284 

247 

139 

992 

394 

45 

230 

9 

285 

231 

140 

969 

395 

61 

231 

4 

286 

264 

341 

973 

396 

60 

232 

7 

287 

301 

142 

949 

397 

43 

233 

3 

288 

315 

143 

967 

398 

52 

234 

8 

289 

307 

344 

963 

399 

26 

235 

0 

290 

373 

345 

919 

400 

45 

236 

21 

291 

384 

346 

922 

401 

34 

237 

11 

292 

359 

347 

845 

402 

13 

238 

23 

293 

395 

348 

866 

403 

13 

239 

22 

294 

424 

349 

823 

404 

16 

240 

22 

295 

424 

350 

775 

405 

55 

241 

20 

296 

481 

351 

719 

406 

0 

242 

20 

297 

468 

352 

721 

407 

4 

243 

14 

298 

512 

353 

728 

408 

22 

244 

16 

299 

506 

354 

701 

409 

16 

245 

12 

300 

526 

355 

672 

410 

5 

246 

26 

301 

581 

356 

672 

411 

30 

247 

27 

302 

580 

357 

638 

412 

0 

248 

20 

303 

622 

358 

571 

413 

1 

249 

13 

304 

643 

359 

563 

414 

12 

250 

24 

305 

720 

360 

595 

415 

0 

251 

38 

306 

702 

361 

515 

416 

15 

252 

24 

307 

713 

362 

527 

417 

0 

253 

24 

308 

787 

363 

469 

418 

4 

254 

24 

309 

816 

364 

458 

419 

0 

255 

22 

310 

799 

365 

408 

420 

0 

256 

28 

311 

818 

366 

398 

421 

0 

257 

31 

312 

851 

367 

398 

422 

0 

258 

48 

313 

847 

368 

324 

423 

9 

259 

37 

314 

866 

369 

364 

424 

7 

260 

68 

315 

919 

370 

364 

425 

7 

261 

57 

316 

897 

371 

344 

426 

16 

262 

40 

1 317 

961 

372 

294 

427 

0 

263 

39 

318 

960 

373 

287 

428 

0 

264 

73 

319 

949 

374 

289 

429 

6 

265 

70 

320 

1036 

375 

259 

430 

0 

266 

83 

321 

976 

376 

256 

431 

8 

267 

86 

322 

1010 

377 

223 

432 

0 

268 

108 

323 

976 

378 

195 

433 

14 

269 

83 

324 

1034 

379 

205 

434 

0 

270 

106 

325 

1085 

380 

172 

435 

13 

271 

92 

326 

1019 

381 

184 

436 

0 

272 

96 

327 

1034 

382 

145 

437 

0 

273 

120 

328 

1050 

383 

133 

438 

0 

274 

112 

329 

1056 

384 

125 

439 

7 

275 

147 

330 

1049 

L 385 

119 

440 

0 



Table I. Continued 


(b) Through fuel (Jet A) 


Channel 

number 

Counts 

Channel 

number 

Counts 

| Channel 
1 number 

Counts 

Channel 

number 


1 

0 

56 

11 

111 

163 

166 

! 19 

2 

0 

57 

3 

112 

172 

167 

38 

3 

0 

58 

5 

113 

165 

168 

42 

4 

0 

59 

10 

114 

168 

169 

40 

5 

0 

60 

5 

115 

162 

170 

53 

6 

0 

61 

0 

116 

195 

171 

27 

7 

0 

62 

4 

117 

176 

172 

32 

8 

0 

63 

1 

118 

169 

173 

43 

9 

0 

64 

3 

119 

178 

174 

38 

10 

0 

65 

0 

120 

163 

175 

35 

11 

0 

66 

8 

121 

175 

176 

30 

12 

0 

67 

10 

122 

155 

177 

38 

13 

0 

68 

2 

123 

157 

178 

29 

14 

0 

69 

6 

124 

171 

179 

51 

15 

0 

70 

4 

125 

182 

180 

26 

16 

0 

71 

0 

126 

173 

181 

37 

17 

0 

72 

0 

127 

164 

182 

32 

18 

0 

73 

4 

128 

181 

183 

30 

19 

0 

74 

0 

129 

140 

184 

36 

20 

0 

75 

6 

130 

155 

185 

34 

21 

0 

76 

4 

131 

151 

186 

29 

22 

0 

77 

3 

132 

153 

187 

23 

23 

0 

78 

5 

133 

128 

188 

29 

24 

0 

79 

8 

134 

147 

189 

33 

25 

0 

80 

6 

135 

122 

190 

24 

26 

0 

81 

0 

136 

141 

191 

10 

27 

0 

82 

6 

137 

121 

192 

33 

28 

0 

83 

9 

138 

113 

193 

19 

29 

0 

84 

2 

139 

no 

194 

22 

30 

0 

85 

4 

140 

93 

195 

8 

31 

0 

86 

4 

141 

99 

196 

20 

32 

0 

87 

1 

142 

89 

197 

14 

33 

0 

88 

3 

143 

104 

198 

27 

34 

0 

89 

8 

144 

81 

199 

7 

35 

0 

90 

10 

145 

74 

200 

12 

36 

0 

91 

0 

146 

87 

201 

14 

37 

0 

92 

0 

147 

55 

202 

6 

38 

0 

93 

0 

148 

71 

203 

15 

39 

0 

94 

0 

149 

54 

204 

22 

40 

0 

95 

5 

150 

57 

205 

0 

41 

0 

96 

0 

151 

33 

206 

8 

42 

0 

97 

0 

152 

22 

207 

11 

43 

0 

98 

0 

153 

46 

208 

7 

44 

2 

99 

13 

154 

30 

209 

15 

45 

0 

100 

0 

155 

36 

210 

18 

46 

5 

101 

134 

156 

35 

211 

16 

47 

2 

102 

140 

157 

25 

212 

9 

48 

1 

103 

114 

158 

29 

213 

3 

49 

3 

104 

142 

159 

33 

214 

6 

50 

0 

105 

160 

160 

22 

215 

4 

51 

7 

106 

133 

161 

39 

216 

10 

52 

8 

107 

142 

162 

31 

217 

4 

53 

9 

108 

149 

163 

54 ! 

218 

2 

54 

0 

109 

186 

164 

32 j 

219 

0 

55 

2 1 

110 

151 

165 

35 

220 

0 



Table I. Concluded 


(b) Concluded 




Table II. Summary of Mass Attenuation Coefficients 
of Common Salt Solutions 


Solution 

number 

Salt solution composition 

p (experimental) , 
cm 2 /g 

^(calculated),* 

cm 2 /g 

1 

100 percent saturated solution 
(35.14 g of salt per 100 cm 3 
of H 2 0) 

0.2242 ± 0.0021 

0.2243 ± 0.0013 

2 

80 percent saturated solution 
(28.11 g of salt per 100 cm 3 
of H 2 0) 

0.2191 ±0.0020 

0.2189 ± 0.0014 

3 

60 percent saturated solution 
(21.08 g of salt per 100 cm 3 
of H 2 0) 

0.2132 ±0.0019 

0.2128 ± 0.0015 

4 

40 percent saturated solution 
(14.06 g of salt per 100 cm 3 
of H 2 0) 

0.2051 ± 0.0019 

0.2061 ± 0.0016 

5 

20 percent saturated solution 
(7.03 g of salt per 100 cm 3 
of H 2 0) 

0.1989 ± 0.0019 

0.1985 ± 0.0019 


* ^(calculated) values were obtained as follows: 


( p ) solution ^1 ( p ) water 4“ 1V 2 ( p ) common salt 


where W\ and W 2 are fractions of the solution by weight. 




Table III. Summary )f Results 


(a) JP-4 fuel 

File name 

Sample I.D 

Run number 

Data collected at 

Fuel type 

Source 

Airline 

Location 

Delivery date 

Peak centroid location (air) 

Peak centroid location (fuel) 

ROI 

Real time, sec 

Live time, sec 

I o 

lx 

Fuel temperature, °C 

Fuel density, g/cm 3 

Linear attenuation coefficient, cm^ 1 . . . 
Mass attenuation coefficient, cm 2 /g .... 


.... JN28026I.CHN 

026 

1 

. . 09:12 on 27 June 88 

JP-4 

NASA 

I.R.D 

Langley 

02-10-88 

Channel number 329.04 
Channel number 328.51 
Channel number 241 to 421 

1800 

1797 

102 744 ±321 

17 436 ±132 

25.0 

0.7520 

.... 0.1392 ± 0.0015 
0.1851 ± 0.0020 


(b) JP-5 fuel 

File name 

Sample I.D 

Run number 

Data collected at 

Fuel type 

Source 

Airline 

Location 

Delivery date 

Peak centroid location (air) 

Peak centroid location (fuel) 

ROI 

Real time, sec 

Live time, sec 

Io 

I* 

Fuel temperature, °C 

Fuel density, g/cm 3 

Linear attenuation coefficient, cm -1 .... 
Mass attenuation coefficient, cm 2 /g 


JN270251.CHN 

025 

1 

. . . 11:04 on 26 June 88 

JP-5 

NASA 

I.R.D. 

Langley 

02-10-88 

Channel number 330.01 
Channel number 329.21 
Channel number 241 to 421 

1800 

1797 

102 679 ±320 

15 756 ±126 

25.0 

0.8086 

0.1492 ± 0.0016 

0.1845 ± 0.0019 
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Table III. Concluded 


(c) Jet A fuel 


File name 

Sample I.D 

Run number 

Data collected at 

Fuel type 

Source 

Airline 

Location 

Delivery date 

Peak centroid location (air) . . 
Peak centroid location (fuel) . . 

ROI 

Real time, sec 

Live time, sec 

lo 

I* 

Fuel temperature, °C 

Fuel density, g/cm 3 

Linear attenuation coefficient, cm -1 
Mass attenuation coefficient, cm 2 /g 


JN270271.CHN 

027 

1 

. . 14:08 on 26 June 88 

Jet A 

NASA 

I.R.D. 

Langley 

02-10-88 

Channel number 330.01 
Channel number 328.80 
Channel number 241 to 421 

1800 

1797 

102 679 ±320 

15 798 ±126 

25.0 

0.8092 

0.1490 ± 0.0016 

0.1841 ± 0.0019 
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Source 

assembly 


Detector 



x glass - 2xL-| = 0.680 cm 
x fluid = l 2 = 10 062 cm 
L3 = 2.540 cm 

x air1 = 2 x L3 = 5.080 cm 
Xg| = L4 = 0.079 cm 

x air2 = x fluid(air) + x air1 = 15.142 cm 
S = americium 241 gamma source 

Figure 1. Geometrical details of fuel cell and associated shields/collimators. 
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Figure 2. Photograph of experimental system. 


L-88-7755 
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100 200 300 400 

Channel number 


Figure 3. Gamma ray spectrum through air. 



Channel number 


Figure 4. Spectrum of gamma rays transmitted through 


Channel 421 


Counts 600 - 



Figure 5. Spectrum of gamma rays transmitted through Jet A fuel. 



Concentration, C (fraction of saturation value) 


Figure 6. Experimental and calculated values of mass a bsorption coefficient of common salt solution 
as a function of concentration. 


in water 



Appendix A 


Flowchart of Program (PGRMAC) 
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Appendix B 

Listing of Computer Program (PGRMAC) 


1 C THIS PR06RAM PRINTS ft SPECTRUM FILE FROM THE EMULATOR 

2 C 

J 1NTE6ER*2 TYPE, MCA, SEE, STRTCH,LNGTDT,SPCOUI (64) 

4 INTEGER*2 SPCINI (64) , BEGREC , ENDREC 

5 INTEGERS SPC I N ( 32 ) , LVETME , RLT I ME , AREA , DAREA , I DSMPL 

6 CHARACTER*! SRTTME 14) ,SRTSEC(2) ,SRTDTE<8> .OUTPUT 130) 

7 CHARACTERS CRCTRL , ANS ,NRUN ,NSttPL*3 

B CHARACTER FUTP *8 1 S0RC»1 0 ,L0CA* 1 0 f DLBT*1 0 , AlfiLN*8 

9 CHARACTER*30 QUT30,F1LNM*I5,FN1*15,FN2*15 

10 C0HH0N/PART1/ SRTTME , SRTSEC , SRTDTE ,FM1 ,FN2,NRUN, NSMPL, MAY! LUID 

11 C0MM0N/PART2/ RLTIME, LVETME, MAXCHN.MAXSPN, LBEG .LEND , MAX All 

12 COMMON/PART3/ ATCFL , ftTCFM ,MCHN (700) ,NCNT (700) ,DIO,DIK,DML DMM 

13 C0MM0N/PART47 ftREA ,DNSTY , 10 , 1 X .POWER , LUNINN, LUNOUT , I REC , THP 

14 C0MN0N/PART5/ FUTP , SORC ,L0CA ,DL0T 

15 EQUIVALENCE (SPCOUT ,SPC0UI 1 ,(SPC1N, SPCINI) 

It EQU1VALENCE10UT30, OUTPUT) 

17 REAL MUA1R.MUGLS ,MUAL, MUFUEL ,MUMAX ,MUM1N,MAXCHN,MAXA1R ,MA FLUID 

18 DATA LUNCON/O/, OUTPUT/30*' 7 

19 DATA CRCTRL/ ‘17, IER/O/ 

20 C 

21 C START 

22 C 

23 LUNINN-1 

24 LUN0UT=3 

25 DPEN117,FILE=’PRN' ,STATUS= ' NEW ’ ) 

2S WRITE (LUNCDN , 90) 

27 90 FORMAT!/,' m**m*mm*m«**m**ft*«**m»m****f »»**' i 

2B | 

29 WRITE (LUNCON, 100) 

30 100 FORMAT(/20X, SPECTRUM PRINT ROUTINE ,//lX, 

31 « ENTER THE DESIRED FILE NAME FROM MCA: 7) 

32 READ (LIJNCON , 5,ERR=B010) FILNM 

33 5 FORMAT IA15) 

34 C - 

35 NSMPL =F I LNM ( 5 : 7 ) 

36 NRUN-F1LNM(8:8) 

37 FN1 (1:4)=FILNM(1:4) 

3B FN1(5:15)='AIR1.CHN' 

39 FN2(1:B)=FN1 (1:8) 

40 FN2<9: 15) = . DAT ‘ 

41 C 

42 CALL CONVERT (NSMPL, IREC) 

43 OPEN ( 19 , F ILE* ' FUELDATA.DAT * ,STATUS= ' OLD' , ACCESS? 'DIRECT ’ 

44 I FORM= ' FORMATTED ' ,RECL=70) 

45 READ ! 19 ,25,REC=IREC) IDSMPL, DNSTY.TMP, FUTP, SORC, LOCA, DLD , 

46 I AIRLN.ADTV 

47 25 FORMAT (I3,lX,F6-4,lX,F4.1,lX f A5,lX,A10,lX,A10,lX,A8,lX,Aj, 

4B * 1X.F5.2) 


49 C - - 

50 OPEN (LUNINN, F ILE*FN1 , STATUS* 'OLD' , ACCESS* 'DIRECT 1 , 

51 I RECL=32) 

52 OPEN ILUNOUT , FILE=FN2, STATUS* NEW) 

53 READ (LUNINN, REC=1 1 TYPE, MCA , SE6 , SRTSEC , RLT IME , LVETME , 

54 < SRTDTE , SRTTME f STRTCH ,LNGTDT 

55 IF (TYPE ,NE. -1) GO TO 8010 

56 C 

57 CALL CALCNTSILUNINN.LUNQUT) 

58 C 

59 CALL CALAREA (LUNOUT ) 

60 MAXA1R=HAXCH« 

61 C 

62 IAIR*AREA 

63 MUAIR=0. 0002132 

64 NUGLS=0. 46654 

65 HUAL*0.6639 

66 XAIR1*5.0S 

67 XAIR2=1 0. 062+X A IRI 

68 X6LS=0.68 

69 XAL=0.079 

70 POKER* IMUAIR*XA!R2)+<MUGLS»XGLS)*iMUAL*XAL) 

71 IO*IAIR*EXP(PONER! 

72 DIO=SQRT (REAL U0>) 

73 C — - 

74 C CHANGE FN1.FN2 FOR EACH SAMPLE 

75 C — 

76 FN1 (1:B1=FJLNH(1:8) 

77 FN1(9:15)=’.CHN' 

78 FM2 ( 1 : 8 ) =F 1 LMM (1:8) 

79 FN2(9: I5>= ' . DAT 

80 LUNINN=11 

81 LUNOUT- 13 

82 OPEN 1LUN1NN, F 1 LE=FN1 , STATUS 1 'OLD , ACCESS 1 'DIRECT' ,RECL=32) 

83 OPEN ILUNOUT ,FILE=FN2, STATUS= NEK ' ) 

84 READ CLUMINN, REC=1 > TYPE, MCA ,SE6,SRTSEC ,RLT IME,LVETME,SRTDTE, 

85 t SRTTME, STRTCH, LNGTDT 

96 IF (TYPE .NE. -1) GO TO 8010 

87 C - — 

88 CALL CALCNTSU.UN1NN, LUNOUT) 

89 C 

90 CALL CALAREA (LUNOUT) 

91 MAXFLU1D=MAXCHN 

92 c - 

93 1 X=AR£A 

94 D 1 X=SORT (REAL (IX)) 

95 RAHAX*R£AL 110*010) /REAL ( i X -0 1 X ) 

96 RANIN=REAL ( 1 0-D10) /REAL TIX+DIX) 

97 PQNERl=MUAlRt(XAlRl)+HUGLS»XGLS*MUAL*XAL 

98 XFUEL= 10. 062 

99 HUMAX= 1 1 . O/XFUEL) * (LOG (RAHAX ) -POKER) 

100 MUHIN* ( 1 . 0/ XFUEL ) * (LOS (RAMIN} -POKER) 
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101 HUHJfct 'IMliHiWHUMlN) '2-0 

10? liH 1 -HUF UEL - HUMAX 

103 DH2*HUFUEL-HUHIN 

104 0ML=SQRT IDMI**2,0+DM2*«2.0) 

105 DMM=DHL/DNSTY 

10b ATCFL=MUFUEL 

107 ATCFH=ATCFL/DNST V 

108 C 

109 CALL PRNTI 

110 C 

111 WRITE (LUNCON, 50) 

112 50 FORMAT!'/, DO YOU WANT TO PRINT OUT CURRENT DATA IY/N)? ’M 

113 RE AD I LUNCON , 55 , ERR=80 1 0 ) ANS 

114 55 FORMAT (All 

115 IF (ANS .Ed. T) THEN 

lib CALL PRNOUT 

117 ELSE 

110 GO TO 999 

119 ENDIF 

120 999 CLOSE (LUN1NN) 

121 1000 STOP 

12? 8000 lER’IER+l 

123 8010 1ER=1ERH 

124 WRITE(LUNC0N,B1K» IER 

125 8110 FORMA) 1 PRINT ERROR* ,14) 

126 GO TO 1000 

127 END 
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Subroutine PRNOUT 


129 

130 
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137 
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166 
167 
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169 


SUBROUTINE PRNOUT 

INTEGERM LVE THE ,RLT IME , AREA , DAREA 
CHARACTER*! SRTTME (4) ,SRTSEC 12) , SRTDTE IB) 

CHARACTER *1 NRUN,NSHPU3 I FN1»I5,FN2»15 
CHARACTER FUJP*B,S0RC*10,LUCA«10,DIDT*10,A!RLN*8 
COMMON/PAR T1 / SRTTME , SR TSEC, SRTDTE , FN1 ,FN2 I NRUN,NSHPL,NAXFLUID 
C0MM0N/PART27 RL T I ME ,L VE THE ,MAXCHN, MAXSPN , LBES, LEND , MAX AI R 
CONNON/PART3/ ftTCFL , ATCFM , NCHN ( 700 ) , NCNT ( 700) , DIO, DI X , 0ML , DMM 
C0MN0N/PART4/ AREA,BNST y , 1 0, 1 X ,P0*4ER , LUNINN, LUNDUT t IREC , TMP 
COMMON/PARTS/ FUTP , SORC ,LOCA , DLDT 
REAL HAXCHN , MAX A1 R , MAXFLUI D 

WRITE 1 17,*> CHAR (14), ' AVIATION FUEL STUDIES 

WRITEt 1 7,10) FN1 

10 FORMAT l///5X,Ht*»mmmHm«#«*»mt*imm»«*mm«« 

« imurntmimim ,//5X, FILE NANE ;',A15) 
NR1TE(17,20) NSMPL,NRUN, SRI THE, SRTSEC, SRTDTE 
20 FORMAT I/5X, SAMPLE 1.0. s * ,A4,//5X, RUN NUMBER :',A2, ,, 

# 5X , DATA COLLECTED AT : ' ,2A1, ' ,2Al , V ,2A1 ,2X, - 

» ' ON ' ,2A1, - ,3A1, -*,3AI) 

HRITE<17,25) FUTP , SORC , LOCA, DLDT 

25 FORMAT I /5X, 'FUEL TYPE : ',A8,//5X, ’SOURCE ; t A10, 

I 7/5X, LOCATION ,A10,//5X, 'DELIVERY DATE , 

# A10) 

WRITE (17,30) MAXAI R, MAXFLUI D , LBEG ,LEND,MAXSPN 
30 FORMAT T/5X, 'PEAK CHANNEL (AIR) ; ‘ ,F7.2,//5X, 

* PEAK CHANNEL (FUEL) s ’ ,F7.2,//5X, R, 0. I. : ,]5, 

* ' ... ',I5,//5X, PEAK COUNTS (FUEL) : , 15) 

WR1 TE ( 1 7, 35) RLT!ME,LVETME,!0,DI0, 1X,D1X 

35 FORMAT (75X, REAL TIME : ' , 15, seconds', / /5X , 

I LIVE TIME i ,15, seconds', 

I 7/5X, 10 : ,18, */-',F5.0,//5X,'IX : ,18, ♦/•• , 

4 F5.0) 

NRITE(17,45) FMP 

45 FORMAT 1/5X, FUEL TEMPERATURE :',F7.1,' [') 

WRI TE ( 1 7 , 40) DNSTy ,ATCFL,DML ,ATCFM,DMM 
40 FORMAT I/5X, FUEL DENSITY : ,F9.5, ' lg/c.3> * ,//5X, 

« LINEAR ATTENU. COEFF. : ,F9.5, ' ♦/-’,F8.5,' (]/«) , 

I //5X, HASS ATTENU. COEFF. : ,F9.5, ' +/-’ 1 F8.5,' (c«2/q) ') 
RETURN " 

END 


V 
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Subroutine PRNT1 


172 BUBROUt 1 NE PRN11 

17J INTEGERS LVETME ,RLT IME , AREA ,DAREA 

174 CHARACTER*! 5RTTME14) .SRTSECI2) ,SRTDTE(B> 

175 CHARACTER*! NRUN,NSMPL*3,FN1*15,FN2»15 

17fc CHARACTER FUTP*B , SORC* 1 0 , L OC A* 10,DLDT« ! 0, A1RLN*B 

177 COMHOM/PART i / SRT THE ,SRTSEC , SRTDTE ,FN1 , FN2 ,NRUN,NSMPL ,HAXFL JID 

178 C0MM0H/PART2/ RLT1HE ,LVETHE ,HAXCHN,MAXSPN,iBE6,LEND,MAXAlR 

179 COMHON/PAR 13/ ATCFL , A fCFM , NCHN < 700) , NCNT (700) ,DI0,DIX,DML,[MH 

180 C0HM0H/PART4/ AREA , DHSTY , 10, IX .POWER , LUNINN , LUNOUT , IREC , THE 

181 C0MM0N/PART5/ FUTP , SQRC ,LOCA , DLDT 

102 REAL MAXCHN,MAXA1R,HAXFLUID 

183 C 

104 WR IT E ( 0 , 0 ) 

185 5 FORMAT I30X, AVIATION FUEL STUDIES /5X ,' ****************** , 

l8i| | ' IHIHHIHHHIHIIIIHIHIIHHIHHHIHIHHHt ‘ ) 

187 WRITE (0,10) FN1 

188 10 FORMAT (5X , 'FILE NAME ;’,A15) 

189 WR ITE (0 , 20> NSHPL ,HRUN,SRTTME ,SRTSEC ,SRTDTE 

190 20 FORMAT 151, SAMPLE I.D. j ’ ,A4,/5X, RUN NUMBER :',A2,/, 

191 « 5X, DATA COLLECTED AT : ' ,2A1, ' S ‘ ,2A1 , ' s ' ,2A1 ,2X, 

192 I ' ON ' ,2A1 , , 3A1 , ,3A1) 

193 WR1TE(0,25> FUTP, SQRC.LOCA, DLDT 

194 25 F0RHAH5X, FUEL TYPE i ' ,A8, /5X, SOURCE i ,A10 

195 i /5X, LOCATION ! ' ,010,751, 'DELIVER! DATE : ,A10) 

196 WR1 TE (0 ,30> MAXAIR , HAXFLU ID ,LBEG ,LEND,MAXSPN 

197 30 FORMAT (SI, PEAK CHANNEL(AIR) :',F7.2,/5X, 

198 I 'PEAK CHANNEL (FUEL) ; ' ,F7 . 2 , 75X , ’ R. 0. I. 

199 I ' ...' ,15,751, 'PEAK COUNTS(FUEL) : ' , 15) 

200 RLTIME-RLTIME/50 

201 LVETME=LVETME/50 

202 WRITE (0,35) RLTIME, LVEIME, 10, DIO, 11,01* 

203 35 FORMAT 151, 'REAL TIME ,15,7 ,51, 'LIVE TIME :',I5,/51, 

204 I 10 i ',18,' 

205 I F5.0,/5X, IX i ',18,' +7- ' ,F5.0) 

206 WRITE LO, 45) TMF 

207 45 FORMAT 151, FUEL TEMPERATURE i ,F7.i> 

200 WRITE 40,40) DNSTV, ATCFL, DML,AFCFM,DMM 

209 40 FORMAT (51, FUEL DENSITY s‘,F9.5,/5X, 

210 « LINEAR ATTENU. COEFF. :',F9.5,' +7 - ' ,F8. 5, 75X , 

211 « MASS ATTENU. COEFF. :,F9.5,’ */-',FB.5> 

212 RETURN 

213 END 
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Subroutine CALCNTS (LIN,LOUT) 


1 
1 
1 

2 
2 
2 
2 
2 
i 


215 

216 

217 

218 

219 

220 
221 
222 
22J 

224 

225 

226 

227 C 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 420 

240 410 

241 500 

242 400 

243 450 

244 8000 

245 


SUBROUTINE CALCNTS < LI N, LOUT) 

INTEGER»2 BEGREC.ENDREC 

1NTEGERM LVETHE,RLT1HE,AREA,DAREA,SPC1N(32) 

CHARACTER # 1 SR T F NE ( 4 ) , SRTSEC ( 2 ) , SR TD TE ( 8 ) 

CHARACTER*! NRUN , NSHPL « 3 , FN1 # 1 5 , FN2* 1 5 
CH AR AC TER FU TP *8 , SORC » 1 0 , LOC A* 1 0 , DLDT » 1 0 , A I RLN»8 
C0HHQN/PART1/ SRTTHE, SRTSEC, SRTDTE,FN1,FN2, NRUN, NSHPL, MAXFLUID 
C0HM0N/PART2/ RLT I HE , LVETHE , HAXCHN , HAXSPN , LBEG , LEND , HAXA IR 
C0HM0N/PART3/ ATCFL , ATCFH, NCHN (700) , NCNT ( 700) , DIO, DIX , DHL , DHH 
C0HH0N/PART4 / AREA , DNSTY , 10, 1 X , PONER , LUN t NN , LUNOUT , I REC , THP 
C0HH0N/PART5/ FUTP,SORC,LOCA,OLOT 
REAL HUAI fi , HUGLS , HUAL , HUFUEL , HAXCHN 

1CHNNL-0 
CHANL I = I CHNNL - 1 
LCHNNL=510 
CHANLL*LCHNNL-1 
BEGREC=CHANL1/B. 

ENDREC=CHANLL/8. 

DO 450 l r BE6REC+2,ENDREC+2 

READ (LIN,REC=I ,ERR=8000) (SPCIN(K) ,K»1 ,8) 

KCHNL=B*( 1-2) 

DO 400 J>1,8 

IFtKCHNL .GT. 1000) GO TO BOOO 
WRI TE (LOUT ,410) KCHNL ,SPCIN ( J ) 

FORHAT (IX, 15, 19) 

KCHNL *KCHNL*1 

CONTINUE 

CONTINUE 

RETURN 

END 
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Subroutine CALAREA (LOUT) 


247 SUBROUTINE CALAREAIL0U1 > 

248 INTEGER*2 BEGREC.ENDREC 

249 INTEGERS LVETHE , RLT 1ME , ftREA , DftREA , SPC1H (32 > 

250 CHARACTER*! SRTTHEI4) .SRTSECI2) ,SRTDTE(8> 

251 CHARACTERS NRUN,NSHPL*3,FN1*15,FN2*15 

252 CHARACTER FUTPt8,SQRC*10,LOCAtlO,OLOT»10,AIRLN#8 

253 COMMON2P ART 1 / SRTTME , SRTSEC , SRI DTE ,FN1 ,FN2 ,HRUN t NSHPL ,t'AXFLUID 

254 C0HHQN/PART2/ RLT IME ,LVETHE , HAXCHN , HAXSPN , LBEB , LEND , HA) A IR 

255 COHH0N/PART3/ ATCFL , ATCFM, MCHBH700) ,NCNT(700) , DIO , DI X ,EML ,OHM 

256 C0HH0N/PART4/ AREA , DNSTY , 10 , 1 X ,POWER f LUNINN ,LUNOUT , IRE( ,THP 

257 C0HH0N/PART5/ FUTP,SORC,LOCA,DLDT 

258 REAL HUA 1 R , MUGLS , HUAL , HUFUEL , HAXCHN , NY 

259 C 

260 REN I NO LOUT 

261 003001=1,500 

1 262 READ (LOUT , 3 10 ,END-300 > NCHNII) ,NCNT(I) 

1 263 310 FORMAT (t X , IS* 19) 

1 264 300 CONTINUE 

265 LBEG=241 

266 LEND=421 

267 AREA=0 

268 HY=0. 

269 BEG1=REAL(NCNT ILBEGS) ) 

270 BEG2=REAL(NCNT ILBEG+2) > 

271 BE63*REAL(NCNf (LBEG*3) ) 

272 HBEG= (BE61 +BES2+BEB3) /3. 0 

273 END1=REAL(NCNT (LEND*!) ) 

274 END2=REAL(NCNT(LEND1) 

275 END3=REAL <NCHT (LEND-11 ) 

276 HEND=lENDl*END2+END3)/3,0 

277 HAV6»0.5*(H8EG*HENDI 

278 DCHN=REAL(LEND-LBE61 

279 SLOPE 1 IHEND-HBEG) /DCHN 

280 00 320 L=LBE6+1 .LENDS 

1 281 DAREA=NCNT(L) 

1 282 AREA=AREA+DAREA 

1 283 H=SLOPE* (REAL (HCHN(L) -LBEG- 1 1 ) *HBEG 

1 284 HY=HY+(REAL(NCNT(L) >-H)»REAL(NCHNIL) I 

1 285 320 CONTINUE 

286 EXCESS*0.5*(HBE6*HEND)»DCHN 

287 MAXCHN*MY/ (REAL(AREA) -EXCESS) 

288 DO 350 I=LBE6+1,LEND*1 

1 289 ERR«REAL (NCHNII) ) -HAXCHN 

1 290 IF(ABSIERR) .LT. 0.51 THEN 

1 291 HAXSPN=NCNTUI 

1 292 GO TO 400 

1 293 ELSE 

1 294 GO TO 350 

1 295 END IF 

1 296 350 CONTINUE 

297 400 RETURN 

298 END 
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Subroutine CONVERT (NSPLJRC) 
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311 C 
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321 


SUBROUTINE CONVERT (NSPLJRC) 

INTEGERM L VETN£,RL TIME, AREA f DAREA 
CHARACTER*! SRTTHE (4) ,SRTSEC(2) ,SRTDTE(8) 

CHARACTER*! NRUN,NSHPL*3 f FNl*15,FN2*15 
CHARACTER FUTP*8,BORC*lO,LOCA*10 1 DLDT*lO t AIRLN*8 
CHARACTER*! CR1 ,CR2,CR3,NSPt*3 

COHHDN/PART!/ SRTTME , SRTSEC , SRTDTE ,FN1 , FN2, NRUN, NSMPL . MAXFLUJD 
C0NN0N/PART2/ RLTIHE f LVETME f HAXCHN t MAXSPN f LBE6,LEND f HAX A IR 
C0HH0N/PART3/ AT CFL t ATCFM , NCHN <700> f NCNT (700) f DIO, DI X . DHL . DHH 
C0MM0N/PART4/ AREA , DNST V , 1 0 , IX f POKER , LUN1 NN , LUNOUT , I REC , THP 
CQMH0N/PART5/ FUTP,SORC,LQCA,DLDT 

CR1 2 NSPL(I: 1) 

CR2=NSPL(2:2) 

CR3=NSPU3:3) 

N1=ICHAR(CR1 ) 

N2*ICHARfCR2) 

N3=ICHAR(CR3) 

iD=(Nl-48}*lOO4(N2“40)*lOf(N3-48> 

IRC 2 ID 
RETURN 
END 
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